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(N ; ABSTRACT 

■ We describe the nieasurement of the beara profiles and window functions for the Atacama Cosmology 

Telescope (ACT), which operated from 2007 to 2010 with kilo-pixel bolometer arrays centered at 148, 
218, and 277 GHz. Maps of Saturn are used to measure the beam shape in each array and for each 
season of observations. Radial profiles are transformed to Fourier space in a way that preserves the 
spatial correlations in the beam uncertainty, to derive window functions relevant for angular power 
spectrum analysis. Several corrections are applied to the resulting beam transforms, including an 
empirical correction measured from the final CMB survey maps to account for the effects of mild 
pointing variation and alignment errors. Observations of Uranus made regularly throughout each 
observing season are used to measure the effects of atmospheric opacity and to monitor deviations in 
telescope focus over the season. Using the WMAP-based calibration of the ACT maps to the CMB 
4^ I blackbody, we obtain precise measurements of the brightness temperatures of the Uranus and Saturn 

Qh. disks at effective frequencies of 149 and 219 GHz. For Uranus we obtain thermodynamic brightness 

^ ■ temperatures T^"^^ = 106.7 ± 2.2 K and T^^^ = 100.1 ± 3.1 K. For Saturn, we model the effects of the 



^ ' ring opacity and emission using a simple model and obtain resulting (unobscured) disk temperatures 

^ ! of Tl"^^ = 137.3 ± 3.2 K and T^^^ = 137.3 ± 4.7 K. 

^,1 Subject headings: cosmology:observations - planets and satellites: individual: Saturn, Uranus 
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1. INTRODUCTION 

The Atacama Cosmology Telescope (ACT) is a 6 m, 
off-axis telescope designed for millimeter wavelength ob- 
servations of the cosmic microwave background (CMB) 
with arcminute resolution (Fowler et al. 2007). The 
telescope is located in Northern Chile, at an eleva- 
tion of 5190 m, where atmospheric conditions are ex- 
cellent for microwave observations from April through 
December. From 2008 to 2010, ACT observed the sky 
with bolometer arrays operating in frequency bands cen- 
tered at 148 GHz, 218 GHz and 277 GHz. The choice 
of bands provides the ability, in principle, to spectrally 
discriminate between signal from the CMB, the Sunyaev- 
Zel'dovich effect (SZE; Sunyaev & Zel'dovich 1970), and 
dusty sources. The 148 GHz array was also in operation 
during a short observing season at the end of 2007. 

The ACT scientific program has produced measure- 
ments of the CMB angular power spectrum at mul- 
tipoles spanning £ ^ 300 to 10000 (Das et al. 2013; 
Dunkley et al. 2013), probing the primary and secondary 
anisotropics (Sievers et al. 2013). ACT has produced 
source catalogs (Marriage et al. 2011a), and galaxy clus- 
ter samples reaching to high redshifts (Marriage et al. 
2011b; Hasselfield et al. 2013). 

The telescope beam acts as a low-pass spatial filter 
on the astrophysical signal from the sky. The effects 
of this filtering are relevant, to varying degrees, to all 
ACT science results. Interpretation of the angular power 
spectrum, in particular, relies on its decomposition into 
components that are correlated across large ranges in an- 
gular scale. The extraction of cosmological parameters 
associated with such components depends on an under- 
standing of the beam shape and the covariant features of 
the window function uncertainty. 

In this work we describe the processing of planet ob- 
servations to obtain the telescope beams, and thus the 
instrument window function, for each detector array for 
each season of operation. We also discuss the use of 
Uranus observations to monitor beam shape and calibra- 
tion. Because the ACT maps are ultimately calibrated to 
the CMB blackbody via the dipole calibration of WMAP 
(Jarosik et al. 2011), we obtain precise measurements of 
the Uranus and Saturn disk temperatures without re- 
sorting to interpolation across unknown source spectra 
or other calibration transfer standards. 

A detailed description of the instrument can be found 
in Swetz et al. (2011), and the main data reduction 
and map-making pipeline is described in Diinner et al. 
(2013). Hincks et al. (2010) present an analysis of the 
ACT beams from the 2008 season; the present approach 
extends this earlier analysis. 

In Section 2 of this paper we discuss the planet 
observations and maps. In Section 3 we obtain the 
telescope window function through an analysis of the 
planet maps and the application of corrections for 
several systematic effects. In Section 4, we interpret 
the apparent brightness of Uranus to calibrate the 
sensitivity of the telescope to atmospheric water vapor 
and focus variation. In Section 5 we use the CMB-based 
calibration of the ACT maps to obtain absolutely cali- 
brated measurements of Uranus and Saturn brightness 
temperatures. 



Table 1 

Number of Saturn and Uranus 
observations selected for calibration 
purposes for each season and array. 





2007 


2008 


2009 


2010 


Saturn 










148 GHz 


22 


22 


13 


35 


218 GHz 




23 


9 


28 


277 GHz 




16 


5 


23 


Uranus 










148 GHz 


16 


37 


94 


113 


218 GHz 




33 


75 


102 


277 GHz 




21 


16 


74 


Note. — 


While 


Uranus was 


visible 



throughout each observing season, Sat- 
urn was available only near the begin- 
ning and end of each season. The number 
of successful planet observations is lower 
at higher frequencies because of increased 
sensitivity in these arrays to atmospheric 
contamination. 

2. OBSERVATIONS AND MAPPING 

2.1. Observations 

The ACT observation strategy is to scan over small 
ranges of azimuth angle while keeping the boresight al- 
titude fixed at 50° to minimize systematic effects due 
to altitude variation. The CMB fields are observed as 
they rise and set through two central azimuth pointings. 
Planet observations are made at this same telescope al- 
titude angle, by briefly interrupting CMB survey scans 
and re-pointing in azimuth to the planet location. 

ACT observed a planet approximately once per night 
of operation in 2007-2010, but only high quality obser- 
vations of Saturn and Uranus are considered in this anal- 
ysis. Saturn and Uranus were preferentially targeted, as 
Saturn's brightness makes it useful for beam profile mea- 
surements, and Uranus' stable brightness is a convenient 
calibration source. While Uranus was visible through- 
out each observing season, Saturn was available only at 
the beginning and end of each season. Jupiter and Mars 
were, respectively, too bright and too low in the sky to 
be useful for calibration work. In Table 1 we list the 
number of high quality Uranus and Saturn observations 
achieved in each season. 

2.2. Planet Maps 

As described in Diinner et al. (2013), maps are made 
from the ACT time stream data according to a maximum 
likelihood technique in which the noise is described in the 
time domain while the signal is assumed to be spatially 
coherent (Tegmark 1997). The bolometer time stream 
data for planet observations are pre-processed in a similar 
fashion to the CMB survey data. The time stream data 
are calibrated, based on detector load curves, after de- 
convolving the effects of the detector time constants and 
low-pass readout filtering. Detectors are then screened 
for quality based on their projections onto a common- 
mode computed using a fiducial set of well-behaved de- 
tectors. Detector screening is performed with the plane- 
tary signal masked out. Maps are made using a dedicated 
code that is optimized for high resolution mapping of sin- 
gle observations; we have confirmed that the main season 
map pipeline code produces compatible results when run 
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on the same data. Each map is solved iteratively, con- 
verging after fewer than 10 iterations. The noise model 
is designed to remove common-mode signal at frequen- 
cies below 0.2 Hz, corresponding to angular scales that 
are much larger than the array size. Because high-pass 
filtering of this kind may suppress signal at large angular 
scales, we study the transfer function of the map making 
solutions in Section 3.4.1. 

2.3. Map Selection 

While planets were observed almost every night, obser- 
vations do not always yield successful maps. The primary 
reason for failed mapping is a low number of live detec- 
tors during periods of high sky loading. This is a greater 
problem at higher frequencies, and for the 277 GHz ar- 
ray in particular leads to many maps with incomplete 
samphng of the planet signal; such maps are discarded. 

Planet observations are sometimes made after sunrise. 
While such observations allow us to study the magnitude 
of pointing and beam focus changes due to thermal de- 
formation of the telescope structure, they are excluded 
from the present analysis. The sunrise cut is also applied 
to CMB observations. 

A small number of maps are cut due to having 
substantially higher noise levels than other maps of 
the same planet for the same array and season. The 
remaining maps are used to establish the telescope beam 
and calibration parameters, as described in the following 
sections. Season mean beam maps, obtained by aligning 
and averaging the selected Saturn maps for each season 
and array, are shown in Figure 1. The Airy pattern 
is easily seen in the maps at 148 GHz and 218 GHz. 
The horizontal streaks are parallel to the telescope scan 
direction. 



3. BEAMS 

In this section we obtain measurements of the tele- 
scope beams for each season and array. Saturn is bright 
and much smaller than the beam solid angle, permitting 
the characterization of our beams with high signal-to- 
noise ratio. Saturn maps are reduced to radial profiles, 
which are then modeled with appropriately chosen basis 
functions and transformed to Fourier space. The result- 
ing transform, the beam modulation transfer function, is 
corrected for a number of systematic effects to produce a 
window function and ^-space covariance suitable for use 
with the ACT CMB survey maps. 

3.1. Radial Profiles 

While the instantaneous ACT beams are slightly (< 
10%) elliptical in cross-section, we analyze them as 
though they were circular, and work exclusively with the 
"symmetrized" (i.e., azimuthally averaged) radial pro- 
files. The resulting beam measurements are suitable for 
any analysis that incorporates a similar simple azimuthal 
average of the signal under study. In particular, the 
cross-linking of CMB survey observations is such that 
the telescope beam contributes to the map, with roughly 
equal weight, at two orientations that differ in rotation 
by approximately 90 degrees; the net effective beam is 
very nearly circular. The goal of our analysis is thus to 
characterize the symmetrized beam. 



Radial profiles are obtained for each Saturn observa- 
tion by averaging map data in radial annuli of width 9'', 
producing average profile measurements i/i at radii 6i. 
On scales smaller than a few arcminutes, the covariant 
noise is sub-dominant to the planetary signal and thus 
the shape of the core beam is very well measured by 
even a single observation of Saturn. However, at larger 
angular scales there is non-negligible contamination from 
the atmosphere. This means that the asymptotic behav- 
ior of the beam cannot be separated from the fluctuating 
background without making some assumptions about the 
beam behavior far from the main lobe. 

The illumination of the ACT optics is controlled by a 
cryogenic Lyot stop at an image of the circular exit pupil; 
this leads to a beam shape described by an Airy pattern 
(for monochromatic radiation) and to the expectation 
that the beam pattern (including the effects of finite de- 
tector size and spectral response) will decay asymptoti- 
cally as 1/^^, where is the angle from the beam peak. 
We fit this model to the radial profile data in order to fix 
the background level of the map and to provide a model 
for extrapolating the beam to large angles. 

The background level of the map and the large-angle 
behavior of the beam are measured by fitting the model 
y{6) = A X (I' /OY + c to the binned profile data in the 
range Oa < Oi < Ob^ with the range chosen such that the 
profile measurements have fallen to below 1% of the peak 
but are still signal dominated. The points i/i with 6i < 
6a, corrected for the background level c and renormalized 
so that y{6 = 0) ^ 1, constitute the "core" of the beam 
profile. The radial profile beyond Oa is referred to as 
the "wing." The fit ranges used for each array, and the 
amplitude of the beam at Oa relative to the peak, are 
given in Table 2. An example of the wing fit for each 
array is shown in Figure 2. 

For each array, the multiple observations of Saturn in 
each season are combined by taking the mean of their 
core profile points and wings. The season mean beam 
profile is thus described by points yi (for each 6i < Oa) 
and a mean wing fit parameter A. The full covariance 
matrix of {^i, ^} is computed from the ensemble 
of Saturn profiles and, as described in the next section, is 
used to propagate covariant error on large angular scales 
into the ^-space beam covariance matrix. 

The solid angle and FWHM of the beam (as well as 
the approximate ellipticity of the non-symmetrized tele- 
scope beam) are shown in Table 2. Solid angles may be 
compared to the values presented in Hincks et al. (2010). 
The reduction in solid angle uncertainty is primarily due 
to our choice to fit the wing over a range of radii closer 
to the beam center. We have doubled the error from its 
formal value to account for systematic variation in solid 
angle results as different fitting ranges are used. 

3.2. Harmonic Transforms and Window Functions 

Sky temperature data may be compared to predic- 
tions based on cosmological models through statistical 
measures such as the angular power spectrum of tem- 
perature fluctuations. The effect of the observing strat- 
egy and telescope beam on the measurement of correla- 
tion functions in CMB data has been discussed in, e.g.. 
White & Srednicki (1995). For ACT, it is important to 
understand the beam and its covariance over a broad 
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Figure 1. Beam maps for each season and array, formed from the average of selected Saturn maps. Gray scale is logarithmic. Each 
panel is 12^ by 12^Each contributing Saturn map incorporates all live detectors in the array. The well-resolved Airy rings indicates that 
the relative detector offsets are accurately known. The frequency labels correspond to the effective frequencies for CMB spectrum radiation 
(Swetz et al. 2011). 



range of angular scales. 

A Gaussian random temperature field T(n) at position 
n on the sphere may be characterized by coefficients Ci 
of the decomposition of its auto-correlation function: 



T(n)T(nO) = ^^C,P,(n-nO, 



(1) 



where are Legendre polynomials. Measurements T(n) 
of the temperature field with a given telescope and ob- 
serving strategy will differ from the true sky tempera- 
ture; for ACT the primary effects are an uncertain cal- 
ibration (which we neglect for now) and the telescope 
beam. For the case of an azimuthally symmetric beam 
that does not vary in time or with position on the sky, 
the auto-correlation function of the observed tempera- 
ture map T(n) will be related to the describing the 
T(n) by (e.g., White & Srednicki 1995) 



(T(n)T(nO 



E 



2^ + 1 
47r 



i 5 



(2) 



where are the multipole coefficients of the beam, de- 
fined by 



{2i+l) 



47r 



P^(cos6'). 



(3) 



In this context, is referred to as the window function. 

The ACT beams are sufficiently concentrated that we 
may work in a flat sky approximation and take ^ to be a 
radial coordinate in the plane. The radial component 
of the Fourier conjugate variable will be £. The az- 
imuthally symmetric two-dimensional Fourier transform 
B{i) of B{0) may be used to obtain the multipole coef- 
ficients according to B^ = B{£) with entirely negligible 
error. We will thus work in the fiat sky approximation, 
but write B£ (with £ in subscript) to emphasize the cor- 
respondence to the multipole coefficients. 

We consider contributions from the wing and from the 
core separately. The transform of the wing (which is 
truncated below Oa) is easily obtained. For the beam 
core, computation of the transform from the binned ra- 
dial profile requires us to interpolate the profile data yi. 
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Table 2 

Summary of beam parameters by array and season. 









Wing 


l^it (93.1) 


Transform Fit (§3.2) 


Beam Properties 




• 
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Array 


Season 


A 


^ D 


B(Oa) 


-^max 


'^mode 


Solid angle 


FWHM 


Ellip. 


Solid angle 










(^lu ; 






(10-9 sr) 


(arcmin) 




(10-9 sr) 


148 GHz 


2007 


A' 


8' 


2.40 ± 0.05 


21900 


11 


224.5 ± 1.4 


1.364' 


1.098 






2008 


A' 


8' 


2.03 ± 0.04 


18400 


9 


216.6 ±1.6 


1.373' 


1.065 


218.2 ±4.0 




2009 


A' 


8' 


2.15 ± 0.09 


16400 


8 


215.8 ±2.6 


1.378' 


1.039 






2010 


A' 


8' 


2.13 ± 0.04 


15800 


9 


215.4 ±1.8 


1.381' 


1.033 




218 GHz 


2008 


2' 


A' 


6.54 ± 0.12 


25500 


6 


116.5 ±1.1 


1.092' 


1.021 


118.2 ±3.0 




2009 


2' 


A' 


7.07 ± 0.29 


25800 


6 


123.9 ±1.7 


1.057' 


1.026 






2010 


2' 


A' 


7.34 ± 0.21 


25100 


6 


125.0 ±1.8 


1.015' 


1.027 




277 GHz 


2008 


2' 


A' 


5.49 ± 0.12 


28400 


7 


98.1 ± 1.9 


0.869' 


1.049 


104.2 ±6.0 




2009 


2' 


A' 


5.22 ± 0.26 


32100 


7 


95.2 ±4.0 


0.879' 


1.118 






2010 


2' 


A' 


5.65 ± 0.18 


28600 


7 


98.9 ±2.0 


0.870' 


1.103 





Note. — Radial profiles are obtained within a distance 6a from the peak; wing and baseline parameters are fit to data 
over the range Oa to Ob- The level B(6a) of the beam, relative to the peak, at the edge of the wing fit region is shown for 
reference. Parameters ^max and nmode describe the basis functions used for fitting a beam model to the radial profile data 
in the beam core. The solid angle and FWHM of the symmetrized instantaneous telescope beam are provided, along with 
the ellipticity (defined as the ratio of major to minor axes of the unsymmetrized mean beam). The solid angles for the 2008 
season may be compared to the values from the independent analysis of Hincks et al. (2010). 
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Figure 2. Wing model fitted to binned radial profile data. Points 
shown, and the best-fit model (solid line) are from a typical indi- 
vidual Saturn observation for each array. A baseline is fit simulta- 
neously and has been removed from the data points. The shaded 
area represents the standard deviation of all individual best-fits to 
wing models for the 2010 season. 



We achieve this by expressing the real space beam as a 
sum of basis functions which we know to have spatial fre- 



quency cut-offs near those determined by the telescope 
optics in each frequency band. 

The primary mirror has a diameter D = 6 m; this 
size limits the optical response to spatial frequencies be- 
low -^max ^ 27rD/X. A natural basis for the Fourier 
space beam is then provided by the Zernike polynomi- 
als R2n{p)i where p = ^/^max- (The Zernike polynomials 
K^{p^(j)) are orthogonal and complete on the unit disk; 
here we impose azimuthal symmetry and need only con- 
sider even n and m = 0.) We thus model the real space 
beam with basis functions f^{0), for non- negative inte- 
ger n, that are proportional to inverse Fourier transforms 
ofijojp): 



no) 



-^2n+l(^^max) 



(4) 



In practice, for some value of £max and a mode count 
nmode the beam model is expressed as 



''^mode 1 



B{e)= Yl 



(5) 



n=0 



where the coefficients are fit parameters. The model 
is fitted to the data by minimizing the of the residu- 
als yi — B{Oi)^ accounting for the full covariance of the yi 
measurements. Values of £max and nmode are adjusted to 
obtain a fit with the minimum number of modes neces- 
sary that gives a per degree of freedom equal to unity. 
The fit yields coefficients and the covariance matrix 
of the errors, {Sa^Sa'^ ). 

Details of the fitting parameters, including ^max, '^mode 
and reduced foi" each season and array, are presented 
in Table 2. The choice of basis functions results in a 
satisfactory expansion of the beam using approximately 
half as many modes as data points. 

After fitting the coefficients in real space, we obtain 
the corresponding beam transform. 



(6) 
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where are the Fourier transforms of f^{0), truncated 
above Oa^ and W£ is the Fourier transform of the extrap- 
olated wing, from 6a to infinity. 

The resulting beam transforms are shown in Figure 3. 
While the beams for the 2008 season are shown, the beam 
features are not substantially different between seasons. 
The figure shows the separate contributions of the core 
and wing to the beam transform. The wing contributes 
significantly only at low i (i.e., below 1000 at 148 GHz, 
and below 2000 at 218 GHz and 277 GHz). 

While this is very similar to the procedure described in 
Hincks et al. (2010), in this work we include covariance 
between radial profile data points, and between the radial 
profile and the wing fit parameter. This results in a 
natural propagation of beam errors on all spatial scales 
into the fitted beam transform and its ^-space covariance 
matrix. Our treatment of beam errors is discussed in the 
next section. 

3.3. Beam and Calibration Covariance 

Using the covariance matrix of the fit coefficients ob- 
tained in Section 3.1, we also obtain the covariance in i 
space: 

+ ^F/'(5a".5w^.)+^(,5u;^<5a")F;}. (7) 

n n 

This covariance includes contributions from the wing fit 
error and the correlation between the wing fit error and 
the coefficients a^. 

This covariance matrix is an £-space representation of 
the covariant features in the radial profiles that con- 
tribute to the B£. However, because the calibration of 
the CMB survey maps is established at a particular an- 
gular scale, we must recast the beam covariance into an 
appropriate form. This is effectively the same procedure 
applied in Page et al. (2003). 

The absolute calibration of the ACT maps at 148 GHz 
is obtained by cross-correlation of the 2008 Southern 
and 2010 Equatorial maps with the WMAP 7- year maps 
(Jarosik et al. 2011) at 94 GHz of the same sky region, 
over angular scales with 300 < £ < 1100 (Hajian et al. 
2011; Das et al. 2013). This leads to an absolute cali- 
bration of the ACT maps centered near i = 700. The 
2009 season Equatorial maps at 148 GHz are then cali- 
brated, over 500 < £ < 2500, to the 2010 season Equa- 
torial maps. The subsequent calibration of the 218 GHz 
maps to the 148 GHz maps is performed in the signal- 
dominated regime with 1000 < £ < 3000. The 277 GHz 
maps have not yet been calibrated to WMAP. For each 
season and array, we factor out the beam amplitude at 
effective calibration scale £ = where L = 700 for 
the 148 GHz array, and L = 1500 for the 218 GHz and 
277 GHz arrays. 

For a celestial temperature signal described by multi- 
pole moments T^,rn, we measure a map M^^rn in detector 
power units, which is related to T by 

Mi^rn = GiTi^rn, (8) 

where Gi = GqBi is the product of the beam (normalized 



such that B£=o = 1) and a global calibration factor Go 
that converts CMB temperature to map units (i.e., to 
detector power units). 

Our calibration to WMAP is a comparison of M and 
T at ^ = L and is thus a measurement of Gl that is 
independent of our beam uncertainty. This leads us to 
recast Equation (8) as 

Mi^m = GLbiT^^m. (9) 

where 

6. = ^. (10) 

Since the errors in and Gl are not correlated, the un- 
certainty in Gi is described by a covariance matrix F^^/ 
that is the sum of a calibration error from the measure- 
ment of Gl, and a term due to correlated error in the 
bf. 

Va' = btbt\^GL? + ^ P,,. - btT.Lv 

-T^tLbf ^btbt'T.LL\^ (11) 

The term in square brackets is the normalized beam co- 
variance that accompanies ACT data releases. The di- 
agonal beam error (i.e., the square root of the diagonal 
entries of the normalized beam covariance) is shown for 
each season and array in Figure 3. At high ^, the error 
in the effective season beam is dominated by an empiri- 
cal map-based correction, which is not included here (see 
Section 3.5). 

While the fitted beam and covariance are an accurate 
description of the binned radial beam profile, they must 
be corrected for a number of systematic effects prior to 
being used to interpret ACT maps. 

3.4. Correction for Systematics 

The beams computed in the preceding section are cor- 
rected for a number of systematic effects that would oth- 
erwise bias the resulting transforms relative to the true 
telescope beam. These are briefly described below. 

3.4.1. Mapping Transfer Function 

Because our map-making procedure includes time do- 
main high-pass filtering, we might expect poor reproduc- 
tion of large spatial scales. We study the mapping trans- 
fer function by injecting a simulated signal into telescope 
time stream data and comparing the output map to the 
input map in Fourier space. The transfer function devi- 
ates significantly from unity only on angular scales larger 
than 20'. The wing fit is performed at somewhat smaller 
radii than this, and we have confirmed that the wing 
fit and extrapolation reproduces the input signal, on all 
scales, to 0.1%. 

A very small (^0.1%) correction due to the 3.5'' pix- 
elization of the planet maps is applied by dividing the 
beam transform by the azimuthal average of the analytic 
pixel window function. 

3.4.2. Radial Binning of Planet Maps 

The binning of the planet map pixel data into annuli 
has a slight impact on the inferred beam transform. This 
is quantified by evaluating the harmonic transform of 
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Figure 3. Upper panels: Beam transforms for each array. Beams shown are for the 2008 season. The total beam is the sum of 
the contribution from the core of the beam, and the extrapolated wing. Lower panels: Diagonal error from the covariance matrix for the 
renormalized beam. The beam normalization has been fixed at calibration t = 700 (1500) for the 148 GHz (218 GHz and 277 GHz) array(s), 
as described in Section 3.3. Note that these curves show uncertainty, not systematic trends, and they do not include additional uncertainty 
from the empirical corrections of Section 3.5. The full covariance matrix shows an anti-correlation between the beam error at angular scales 
above the calibration i and beam error below the calibration i. The window function is and thus its fractional error is 25Bi/Bi. 



data points taken from a model of the radial beam pro- 
file, and comparing it to the harmonic transform of points 
that include simulated binning of the radial beam pro- 
file. This resulting correction is approximately 1% at 
I = 10000. 

3.4.3. Saturn Disk and Ring Shape 

The Saturn angular diameter of approximately 8^' is 
large enough that we do not treat it as a point source. 
For each season and array, we deconvolve Saturn's shape 
assuming that it is a disk with solid angle equal to the 
mean solid angle for all Saturn observations that con- 
tribute to the mean instantaneous beam measurement. 
This correction is approximately 2% at £ = 10000. 

While Saturn's rings complicate the spatial distribu- 
tion of the planetary signal, these effects are negligible 
at ^ < 20000, aside from an overall change in average 
brightness, and need not be considered in the deconvo- 
lution procedure. 

3.5. Mean Instantaneous vs. Effective Beam 

After applying the corrections described in Section 3.4, 
the resulting beam describes the telescope response to a 
point-like radiation source with approximately Rayleigh- 
Jeans spectrum, averaged over the focal plane of each 
array. We refer to this beam as the "mean instanta- 
neous beam," to emphasize that the combination of ob- 
servations taken at many different times may lead to an 



effective window function that is different from the in- 
stantaneous one. 

For use in particular contexts, we compute an effective 
beam that includes corrections for alternative frequency 
spectra, and for pointing variation or other cumulative 
effects resulting from the combination of observations 
taken at many different times. 

We apply a simple first-order correction to obtain the 
effective beam for different radiation spectra. For radia- 
tion with a band effective frequency the beam is taken 
to be 



(12) 



where urj is the effective frequency for radiation with 
a Rayleigh- Jeans (RJ) spectrum. Effective frequencies 
for various spectral types are presented in Swetz et al. 
(2011). 

In the absence of planetary sources, the absolute point- 
ing registration of individual ACT observations is not 
known. Since each pixel of the ACT survey maps con- 
tains contributions from data acquired on many different 
nights throughout each season, the resulting season effec- 
tive beam for these maps is less sharp than the instan- 
taneous beam obtained from carefully aligned individual 
planet observations. 

Pointing repeatability can be estimated from the vari- 
ation in apparent planet positions, relative to expecta- 
tions, in each season and array. As summarized in Ta- 
ble 3, these indicate that the repeatability of telescope 
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pointing is at the cTpianet = 3'' to 8'' (RMS) level. Point- 
ing variation in CMB observations may in principle be 
smaller than cTpianet, since CMB observations are per- 
formed at the same azimuth angles each night, while 
planets observations span a wider range of azimuths. 

In addition to pointing repeatability, the season effec- 
tive beam may also be diluted due to errors in the global 
pointing model, which is used to combine observations 
made at different azimuth angles. These global adjust- 
ments are estimated using point source positions mea- 
sured in (non-cross-linked) maps made from only rising 
or only setting data. 

While pointing variance and global alignment may con- 
tribute to dilution of the season beam, we must also ac- 
knowledge that the beam may change slightly over the 
course of the season in a way that is not captured by 
the Saturn observations. (Changes in the beam over the 
course of the season, and the resulting impact on cali- 
bration, are assessed using Uranus observations in Sec- 
tion 4.) We thus seek an empirical measure of the differ- 
ence between the season effective beam and the instan- 
taneous beam obtained from Saturn observations. 

We parametrize the difference between the season ef- 
fective beam B^^ and the instantaneous beam B£ as a 
Gaussian in i: 

Bf = Bex e-^(^+i)^/2. (13) 

If the correction is interpreted as arising purely from 
Gaussian pointing error, then V is simply the pointing 
variance in square radians. While this motivates the form 
of the correction, our broader interpretation permits V 
to differ from the value expected based on pointing vari- 
ance measured from planets, or to be less than zero if 
the season effective beam is somewhat sharper than the 
ensemble of Saturn observations indicates. 

We measure V for each season and array using bright 
point sources in the season CMB survey maps. Point 
sources with signal to noise ratios greater than 10 are 
identified using a matched filter. A beam model is fit to 
the map in the vicinity of each source, producing a value 
of V. The weighted mean and error of these individual 
fits, for each season and array, are used to form the ef- 
fective beam for use with the survey maps. While this 
point source analysis gives us a secondary probe of the 
beam core, we note that the Saturn observations are still 
key for measuring the behavior in the wings of the beam. 

The values of the season correction parameter V are 
presented in Table 3. The measurement of V is some- 
what susceptible to changes in the fitting parameters and 
thus the quoted errors have been inflated from the for- 
mal values by factors of 2 and 6 for the 148 GHz and 
218 GHz arrays, respectively. While values of V for the 
148 GHz array are very similar between seasons, values 
for 218 GHz are more variable. The analysis of Uranus 
observations in Section 4.1 shows that in 2010 the tele- 
scope beam was more sharply focused, relative to the 
mean beam estimated from Saturn, during most of the 
observing season (see Section 4.1); this results in a neg- 
ative value for V at 218 GHz. 

4. PLANET BRIGHTNESS MEASUREMENTS 

In this section we use observations of Uranus and Sat- 
urn to study the variation in the system gain due to 
changes in water vapor level and telescope focus. 



Table 3 

Pointing variance and beam correction parameters. 



Season 


Array 


planet 

(arcsec^) 


V 

(arcsec^) 


Attenuation 
at i = 5000 


2007 


148 GHz 


53 ± 12 






2008 


148 GHz 
218 GHz 
277 GHz 


24 ± 4 
23 ± 5 
11 ± 2 


25± 10 
143 ± 51 


0.993 ± 0.003 
0.962 ±0.013 


2009 


148 GHz 
218 GHz 
277 GHz 


58 ± 9 
43 ± 8 
49 ±18 


26 ± 6 
46 ±27 


0.993 ± 0.002 
0.988 ± 0.007 


2010 


148 GHz 
218 GHz 
277 GHz 


31 ± 3 
18 ± 2 
25 ± 3 


29 ± 6 
-5 ±23 


0.992 ± 0.002 
1.002 ± 0.006 



Note. — Pointing variance ^^p^^j^et measured from 
planet observations, and the season effective beam correc- 
tion parameter V is measured from point sources in the full- 
season CMB survey maps, where available. The parameter 
V is used to correct the high--^ beam for use with the survey 
maps; the resulting attenuation or inflation of the beam at 
£ = 5000 is provided for reference. The covariant error due 
to this correction is included in the season effective beam 
covariance matrix. 

4.1. Planet Amplitudes 

We characterize the apparent brightness of a planet in 
a map by examining the ratio of the Fourier transform of 
the map to expectations based on the beam and planet 
shape. While ultimately equivalent to procedures that 
involve measuring source peak heights in filtered real- 
space maps, the Fourier space approach emphasizes the 
role of spatial filtering in providing precise planet ampli- 
tude measurements. It also allows us to quantify differ- 
ences between the planet shape and expectations based 
on the beam measurements in a physically meaningful 
way. 

Once again working in the flat sky limit, we let 
represent a vector in the plane of the sky centered on 
the planet. The Fourier conjugate variable is denoted 
as £ and is written in subscript. The signal from a 
planet is modeled as a circular disk of uniform RJ tem- 
perature . For each planet observation, we use the 
total disk solid angle to compute angular radius 
6^ = (I^^/tt)^/^. The Fourier transform of the plane- 
tary signal has components 

Te = T^^^f , (14) 

where nf = n^x {2Ji{£0^)/{£0^)) is the Fourier trans- 
form of a disk of radius 6^ . For Uranus, ^ at 
the ACT angular scales; for Saturn, the finite disk size 
cannot be neglected, but the effects due to its oblateness 
are negligible. Planetary radii are taken at the 1 bar at- 
mospheric surface, as reported in Archinal et al. (2011). 
The Saturn equatorial radius is 60268 ± 4 km and the 
polar radius is 54364 ± 10 km; the Uranus equatorial ra- 
dius is 25559 ± 4 km and the polar radius is 24973 ± 20 
km. 

Adapting equation (8) to the flat sky, the ACT map 
M{0) of the planetary signal will have Fourier compo- 
nents 

M^ = GoBiT^n^. (15) 
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For each planet map, we compute the amplitude 
of the planet at multiple bins in I. The map is centered, 
apodized, and Fourier transformed. To suppress atmo- 
spheric contamination, which contributes faint horizontal 
streaks to the map, we mask Fourier components having 
l^a^l < 1700^. The remaining complex Fourier amplitudes 
are averaged inside annuli of width 600 and the 

imaginary part is discarded. We then compute the ratio 
of Mi to BiQ.^ (which has been corrected to account for 
the apodization and masking applied to obtain M^). 

The result is a curve of planet amplitude measure- 
ments. 



148 GHz - 2008 



218 GHz - 2010 



(16) 



for 2000 < ^ < 14000. In the absence of any gain or 
beam variation, should be equal to GqT^ ^ the appar- 
ent brightness of the planet in detector power units. 

The mean and variance of ai curves for selected seasons 
and arrays are shown in Figure 4. The ai curves are not 
in all cases consistent with a constant value. This may be 
attributed to variations in the focus from night to night. 
To parametrize these deviations and to make possible an 
assessment of the impact they might have on the system 

calibration, we fit parameters a and P for each curve 
according to the following model: 



a 



5000/ ' 



(17) 



The mean angular scale I is taken as the center of the fit- 
ting range. The parameters a and p represent the mean 
amplitude of the planet and the deviation from expected 
focus of the beam, respectively. For example, P = 0.01 
would indicate that the the planetary signal is stronger 
by 1% at £ 10000 relative to I ^ 5000, given the mean 
season beam, and thus the telescope beam at the time 
of observation was slightly more sharply focused. We 

would expect observations associated with negative j3 to 
also have a smaller a than expected, since a less sharply 
focused beam will have poorer overall efficiency. 

The extraction of amplitudes a and focus parameters 
P allows us to quantify deviations of focus over the sea- 
son, and to estimate the effects of such deviations on the 
overall system calibration. This is accomplished in the 
next section, where a model for system calibration is ob- 
tained from the brightness and focus measurements of 
Uranus. 

4.2. Calibration Parameters from Uranus Observations 

The measurements of a and P from Uranus observa- 
tions, along with measurements of the precipitable water 
vapor (PWV) level for most of these observations, per- 
mit us to characterize the impact of PWV level and focus 
variation on the telescope calibration. This analysis does 
not require that the brightness temperature of Uranus 
be precisely known, but assumes that the brightness is 
roughly constant over the observing season. (We do not 

^ No such masking or filtering is performed on the Saturn maps 
prior to their use in the determination of the telescope beam; the 
masking is more important for Uranus maps, because the signal is 
weaker. 




4000 



8000 12000 



4000 



8000 



12000 



Figure 4. Variation of planet amplitude measurements with an- 
gular scale. Each panel shows mean fractional deviation of 
relative to a. (solid line) and the standard deviation of all indi- 
vidual OLgloL curves (grey band) for a particular season, array, and 
planet. See Section 4.1 for definitions of otn and S. We show results 
from 148 GHz 2008 (left panels), where there is a high degree of 
consistency between the planet observations and the season mean 
beam. We also show results for 218 GHz 2010 (right panels), where 
there is somewhat more variance in the curves (which is typical of 
the 218 GHz array) and where the Uranus observations deviate 
significantly from the mean beam (which is seen in all arrays in 
2010). Differences between Uranus and Saturn are indicative of 
2% changes in the beam focus during the season; such changes 
have been accounted for in the season effective beam by fitting a 
correction parameter to maps of point sources in the season maps 
as described in Section 3.5. 



perform the same analysis on Saturn, because the effec- 
tive brightness varies with the ring opening angle; see 
Section 5.2). 

The Atmospheric Transmission at Microwaves model 
(ATM; Pardo et al. 2001) provides estimates of atmo- 
spheric emission and absorption at millimeter and sub- 
millimeter wavelengths. The ACT frequency bands avoid 
strong molecular resonances and are thus primarily sus- 
ceptible to continuum emission and absorption due to 
water vapor. For the purposes of calibration we are in- 
terested in the atmospheric transmission in a given fre- 
quency band. The transmission at the ACT observation 
altitude is written as 



exp (-(tw^ + Td)/ sin6>o) , 



(18) 



where w is the precipitable water vapor column density 
at zenith and r^w and Td are the "wet" and "dry" op- 
tical depths at zenith. The parameters and Td vary 
negligibly over the range of temperatures and pressures 
experienced at the telescope site. 

Because planet observations discussed in this work are 
obtained at the same fiducial altitude used for CMB ob- 
servations, the effects of Td are common to both observa- 
tions and cannot be measured from the planet observa- 
tions. We define the overall system gain Go such that it 
apphes at the typical PWV level of = 0.44 mm. The 
i(;-dependent gain factor is then 

G^(w) = Goexp(-(rw t^{w - wq)/ sin6>) . (19) 

Here Tw is as provided by the ATM model, while Tx is 
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a fit parameter that could indicate a deviation from the 
ATM model but is more likely due to a bias in the detec- 
tor calibration procedure that is sensitive to sky loading. 
Estimates of Tw for the ACT bands were computed us- 
ing the ATM model by Marriage (2006); these values are 
provided for reference in Table 4. 

To assess the impact of small variations in telescope 
focus on planet calibration, we include a dependence on 
the difference between the focus parameter /3 (defined 
in Equation (17)) and the season mean focus parameter 
^. Because a planet was observed approximately once 
per night, we take /3 to be the average of measured /3 
for all planet observations for the array and season un- 
der consideration. We adopt an exponential form out 
of convenience; calibration variation is small so a linear 
correction would be equivalent. 

The full model for the Uranus brightness data S, (3 is 
then given by 

exp(-(rw +rx)(^ - wq)/ sin6>) , (20) 

where the fit parameters T^Gq, m, and Tx represent over- 
all calibration, sensitivity to variations in telescope fo- 
cus, and sensitivity to PWV beyond predictions based 
on ATM, respectively. 

This model is fitted to the Uranus brightness amplitude 
data for each season and frequency array. The parame- 
ter values are presented in Table 4. The residuals have 
scatter of approximately 2%, 2% to 6%, and 6% in the 
148 GHz, 218 GHz and 277 GHz arrays respectively; this 
scatter is not explained by uncertainties in planet ampli- 
tude measurement. Values obtained for Tx are in good 
agreement between seasons for the 148 GHz and 218 GHz 
arrays. They are inconsistent with zero (at the ~ 2-<j to 
4-cr level) , which we attribute to systematic detector cal- 
ibration error. For the 277 GHz array there is some slight 
disagreement between the 2009 and 2010 results, which 
we attribute to the low number of observations in 2009. 
The values of Tx are used in combination with to cor- 
rect the time-ordered data for opacity prior to creation 
of full-season survey maps, as described in Diinner et al. 
(2013). 

While for some seasons and arrays we obtain better fits 
to the data by including the parameter m, the impact of 
this additional parameter on inferred overall system cal- 
ibration is small. To show this, we refit the model of 
equation (20) with parameter m fixed to 0. The result- 
ing calibration factor T^Gq^~^^ may be compared to 
the case where m was free to vary. The differences are 
presented in Table 4; the change in overall calibration 
is 0.6% or smaller, and always less than the l-cr error 
on T^Go, with the exception of the 2009 data for the 
218 GHz array. This exceptional case also has a large 
(9%) uncertainty in its absolute WMAP-based calibra- 
tion. 

The fit values of T^Go may be used to provide an ab- 
solute calibration of each array (should be known) 
or to provide a measurement of the Uranus disk temper- 
ature based on an independent calibration of Gq. The 
latter possibility is discussed in the next section, leading 
to the measurements of Uranus and Saturn brightness 



Table 4 

Properties of planetary data and calibration model fit 
parameters by season and array. 



Quantity 


148 GHz 


Array 
218 GHz 


277 GHz 


Tw (mm~^) 


0.019 


0.044 


0.075 


Tx (mm~^) 
2007 
2008 
2009 
2010 


0.016 ±0.010 
0.014 ±0.004 
0.009 ±0.005 
0.012 ±0.003 


0.040 ± 0.008 
0.021 ±0.017 
0.045 ± 0.008 


0.033 ± 0.063 
-0.008 ± 0.041 
0.083 ±0.016 


Error in Gq 
2007 
2008 
2009 
2010 


(%) 

0.9 
0.2 
0.3 
0.2 


0.4 
1.1 
0.4 


1.3 
2.2 
0.8 


m 

2007 
2008 
2009 
2010 


0.7 ±0.6 
0.5 ±0.3 
0.5 ±0.3 
0.8 ± 0.1 


-0.5 ±0.6 
3.4 ±0.5 
0.5 ± 0.2 


1.2 ± 1.3 
1.2 ±0.7 
1.5 ± 0.4 


(3 mean (rms) 
2007 
2008 
2009 
2010 


-0.017 (0.010) 
0.002 (0.006) 
0.004 (0.011) 
0.015 (0.015) 


0.011 (0.010) 
0.008 (0.016) 
0.021 (0.019) 


0.011 (0.015) 
0.007 (0.020) 
0.017 (0.019) 


2007 
2008 
2009 
2010 


- l) (%) 
-0.6 
0.1 
0.2 
0.3 


-0.3 
1.4 
0.2 


0.5 
0.4 

0.4 


Note. — The opacity parameters Tw from the ATM model are 



provided for reference. Fit parameters Tx and m parametrize sen- 
sitivity of Uranus apparent brightness to PWV and focus param- 
eter, respectively. The mean and standard deviation of the focus 
parameter (3 for all planet observations in each season is provided 

for reference. The overall calibration shift Gq"^~^V^o — 1 is the 
systematic calibration change if m is fixed to zero instead of being 
free to vary. 

temperatures presented in Section 5. 

4.3. Relation to WMAP-Based Calibration 

The absolute calibration of the ACT full-season 
CMB survey maps is obtained through cross-correlation 
in -^-space (Hajian et al. 2011) to the 94 GHz maps 
from WMAP7 (Jarosik et al. 2011). As described by 
Das et al. (2013), the WMAP maps are used to calibrate 
the ACT 2008 Southern and 2010 Equatorial field maps 
at 148 GHz directly; the latter maps are then used to cal- 
ibrate the 148 GHz map for 2009 and the 218 GHz maps 
for all seasons. 

While the data entering the survey maps are corrected 
for atmospheric opacity, they have not been corrected 
for any kind of time variation in the telescope beam. In- 
stead, we have obtained an effective season beam suitable 
for use with the survey maps by comparing our mean in- 
stantaneous beam to point sources in the survey maps, 
as described in Section 3.5. This correction has negligi- 
ble effect at low £ where the calibration to WMAP takes 
place. Furthermore, the WMAP calibration should be 
associated with the telescope gain at the season- average 
focus parameter. Thus the calibration factor obtained 
from the WMAP calibration is compatible with Go as 
defined by equation (20). 
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5. BRIGHTNESS TEMPERATURES OF URANUS AND 
SATURN 

Using the WMAP-based calibration of the ACT sur- 
vey maps at 148 GHz and 218 GHz, we convert the Sat- 
urn and Uranus amplitude measurements to calibrated 
brightness temperatures. A geometrically detailed, two- 
parameter model for Saturn disk and ring temperatures 
is fit to the Saturn brightness data. Because CMB survey 
maps for the 277 GHz array have not yet been calibrated 
to WMAP, we do not consider that array in this analysis. 

In what follows, we refer to both the Ray leigh- Jeans 
(RJ) and brightness temperatures of the planet.^ The 
RJ temperature Trj is related to the brightness tem- 
perature Tb by equating the specific intensity of an RJ 
spectrum and a blackbody spectrum at the band effective 
frequency i^eff- 

B,^,,{T^)=2vl^kBT^,/c\ (21) 

where Bjj(T) is the blackbody spectral radiance. The 
spectra of Saturn and Uranus are each sufficiently close to 
the RJ limit that we take the source effective frequencies 
to be the RJ band centers of 149.0 GHz and 219.1 GHz 
(Swetz et al. 2011). 

5.1. Uranus 

Based on our fitted values of T^Gq (Section 4.2) and 
the absolute calibration of Go using WMAP data, we ob- 
tain measurements of T^. Because the calibration factor 
is obtained from the primary anisotropics of the CMB, 
we must account for the different frequency spectrum of 
the planets relative to the CMB blackbody. The effec- 
tive spectral index of the flux of the CMB in the ACT 
bands is 1.0 (0.0) in the 148 GHz (218 GHz) bands, while 
the spectral index of Uranus is approximately 1.65 in 
both bands (based on, e.g.. Griffin & Orton 1993). The 
3.5 GHz uncertainty in the ACT band centers and the 
effective frequencies provided in Swetz et al. (2011) re- 
sult in an additional 1.6% (2.7%) error in the Uranus 
brightness determinations at 148 GHz (218 GHz). 

The Rayleigh- Jeans temperatures obtained for each 
season and array are presented in Table 5. The value 
quoted for the 148 GHz (218 GHz) array is valid at effec- 
tive frequency of 149.0 (219.1) GHz. All RJ temperatures 
have been inflated by 0.5 (0.3) K to compensate for the 
background provided by the CMB in our bands. Bright- 
ness temperatures may be computed from the RJ values 
by adding 3.5 (5.1) K. 

The results are consistent between seasons (even af- 
ter removing the correlated error due to the frequency 
spectrum correction), with uncertainty dominated by the 
overall instrument calibration uncertainty rather than 
statistical uncertainty from planet brightness measure- 
ments. 

From the RJ temperature measurements in each sea- 
son, we compute a weighted mean disk temperature for 
Uranus. The WMAP-based calibration error is corre- 
lated between the 2009 and 2010 seasons, because the 
2009 map is calibrated to the 2010 ACT map, and the 
ACT band center uncertainty is correlated across all sea- 

^ Our planet amplitude measurements include a deconvolution 
of the telescope beam and thus we do not work with the planet's 
"antenna temperature" directly. 



sons. We account for this in the mean and error pre- 
sented. At 148 GHz, we obtain a RJ (brightness) temper- 
ature of 103.2 K (106.7 K), with 2.1% error. At 218 GHz, 
we obtain temperature 95.0 K (100.1 K), with 3.2% error. 

We compare our Uranus temperature measurements 
to the widely used model of Griffin & Orton (1993, 
hereafter G&O). The G&O model provides brightness 
temperatures for Uranus and Neptune in the ^ 100 
to 1000 GHz range based on precise measurements of 
brightness ratios of each planet to Mars. The absolute 
calibration is obtained by interpolating, in the logarithm 
of the frequency, between the Wright (1976) model for 
Mars temperature at 3.5/im (857 GHz), and the (Ulich 
1981) model for Mars temperature at 90 GHz. They 
adopt a 5% systematic error for this Ulich- Wright hy- 
brid model. 

The G&O model predicts a brightness temperature 
of 112 K at 148 GHz, approximately 5% larger than 
our measured value. This suggests that the Ulich- 
Wright model for Mars brightness is 5% high at 148 GHz. 
Weiland et al. (2011) have observed a similar 5% discrep- 
ancy between Mars temperatures and the Ulich model at 
94 GHz. At 218 GHz, the G&O model gives a brightness 
temperature of 99 K, which is 1% smaller than our mea- 
sured value; this difference is smaller than our calibration 
uncertainty. 

Long term studies of Uranus show a variation in mean 
disk brightness at 8.6 GHz (Klein & Hofstadter 2006) 
and 90 GHz (Kramer et al. 2008) on decade time scales. 
If this variation is attributed to differences in the Uranus 
surface brightness at different planet latitudes, it implies 
that the mean disk temperature is roughly 10% larger 
when the South pole, rather than the equator, is directed 
to towards the observer.^ The ACT measurements over 
the 2008-2010 period span Uranus sub-Earth latitudes 
from 1° to 13°. In comparison, the observations pre- 
sented by G&O were made in 1990-1992 and span sub- 
Earth latitudes from -70° to -60°. It is then possible 
that, rather than a difference in calibration, Uranus was 
brighter in these bands at the time the G&O data were 
taken. However, if one accepts the 5% downward recal- 
ibration of the Ulich Mars model at 94 GHz, then one 
should accept a downward recalibration (of 4-5%) of the 
G&O data around 148 GHz. This brings the ACT and 
G&O measurements into agreement, with no evidence for 
any variation in Uranus brightness temperature due to 
sub-Earth latitude. 

Our results at 148 GHz are consistent with the analysis 
of Sayers et al. (2012), who use WMAP measurements of 
Uranus and Neptune at 94 GHz (Weiland et al. 2011) to 
recalibrate the G&O model at 143 GHz. It is important 
to note that while our measurements rely on a detailed 
understanding of the telescope beam, our absolute cali- 
bration standard is the CMB, for which the brightness as 
a function of frequency is known. Our errors thus do not 
include any significant systematic contribution due to in- 
terpolation of unknown source spectra into our frequency 
bands. 

In Figure 5 we show our brightness temperature mea- 
surements along with WMAP measurements below 100 

^ The center of the projected planet disk, as observed from 
Earth, is called the sub-Earth point. The Uranian latitude of this 
point is called the sub-Earth latitude. 
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Table 5 

RJ temperature measurements of Uranus and Saturn. 





T^J (K) 






zio \jYVL 


Uranus (§5.1) 






2008 


101.5 ± 2.0 


93.2 ±2.3 


2009 


103.8 ±2.2 


91.5 ±8.6 


2010 


105.7 ±2.1 


97.6 ±2.6 


Combined 


103.2 ±2.2 


95.0 ±3.1 


Saturn, Wll model (§5.2) 






Tdisk 


133.8 ±3.2 


132.2 ±4.7 


Tring 


17.7 ±2.2 


12.3 ±4.0 


be 1 d.o.f.] 


[7.6 / 5] 


[6.5 / 5] 


Saturn, single T model (§5.2) 






^disk 


131.2 ±3.1 


127.5 ±4.4 


bd 1 d.o.f.] 


[12.0 / 6] 


[11.2 / 6] 



Note. — Values for Uranus (Section 5.1) are shown 
for each season, as well as a weighted mean that takes full 
account of calibration covariance between seasons. Saturn 
results (Section 5.2) are provided for both a two-component 
disk±ring model and a disk-only model. All temperatures 
are RJ and have been corrected for the CMB to indicate the 
brightness that would be seen relative to an empty (T=0) 
sky. To obtain thermodynamic brightness temperatures at 
148 GHz (218 GHz), add 3.5 K (5.1 K) to the RJ values. 

GHz from Weiland et al. (2011) and with data used by 
G&O to fit their model in the 90 to 1000 GHz range. 

As discussed above, systematic error in the G&O 
model is dominated by uncertainty in the Mars spectrum 
below 1000 GHz. The ACT and WMAP data provide 
new, absolutely calibrated temperature measurements 
between 20 and 220 GHz. We thus fit a simple empiri- 
cal model for the temperature using the ACT, WMAP, 
and the higher frequency G&O data. To avoid the strong 
spectral feature at 30 GHz, we include only the two high- 
est frequency bands of WMAP. We use the three data 
points in G&O above 600 GHz. We model the Uranus 
brightness temperature as a function of = (z^/lOO GHz) 

by 



Fk 



-— — = ao + ai log 



10* 



a2 log 



10 ' 



(22) 



and obtain best fit coefficients (ao,ai,a2) = 
(121,-78.2,18.2). 

The ACT, WMAP, and G&O brightness temperature 
measurements are shown in Figure 5, along with the 
G&O model and our empirical model. 

5.2. Saturn 

Measurements of the Saturn disk temperature are com- 
plicated by the presence of the rings, whose dust both 
obscures the main disk and radiates at a lower tempera- 
ture. The ACT observations of Saturn span ring opening 
angles from -2 to 12 degrees and provide an opportunity 
to explore the separate contribution of these two compo- 
nents to the total effective brightness of the planet. 

We follow closely the approach of Weiland et al. (2011, 
hereafter Wll), who apply a two-parameter model in 
which the opacity and geometry of the rings are fixed, 
while the disk temperature and an effective temperature 
for the rings are free to vary. The ring dimensions and 
normal opacities are taken from Dunn et al. (2002). For 
each of the 148 GHz and 218 GHz arrays, we fit Tdisk 
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Figure 5. Measurements of the Uranus brightness tem- 
perature from this work, Weiland et al. (2011, WMAP), and 
Griffin & Orton (1993, G&O). The G&O points include data from 
that work as well as Ulich (1981) and Orton et al. (1986). The 
dotted line shows the model of G&O. G&O data are calibrated to 
a Ulich- Wright hybrid model for Mars brightness that interpolates 
between 90 and 857 GHz and carries an estimated 5% systematic 
error. The solid line is our best-fit empirical model, using only the 
two highest frequency points from WMAP, the two ACT points, 
and the G&O points above 600 GHz. 



and Tring, which describe the mean disk temperature and 
mean effective brightness temperature of the rings (see 
Wll for the detailed definition), respectively. 

The measurements of a and j3 for each Saturn observa- 
tion are converted to RJ temperatures, including correc- 
tions for PWV level and focus parameter, and with the 
CMB-based calibration applied. Because the scatter in 
the temperatures is not fully accounted for by the errors 
in the amplitude measurements and correction parame- 
ters, we combine observations made in single 15 to 30 day 
periods, taking the uncertainty in each combined mea- 
surement to be the error in the mean of the contributing 
data. As a result, a small number of temporally isolated 
observations are discarded. The binned data points, and 
the best fit model, are shown in Figure 6. 

The fit of the Wll model includes a full accounting of 
the non-trivial correlations in the calibration uncertainty 
of the data points. The resulting models are good fits, 
with P^^ degree of freedom of 7.6/5 (6.5/5) in the 
148 GHz (218 GHz) array. The model fit parameters and 
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uncertainties are presented in Table 5. 

For comparison, we also fit a single temperature model 
to the binned Saturn data, finding poorer fits and mean 
temperatures that are 1-cr less than the disk tempera- 
tures obtained in the full model. These temperatures 
and statistics are also presented in Table 5. 
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Figure 7. Measurements of the Saturn disk brightness temper- 
ature from ACT, Weiland et al. (2011, WMAP), and Goldin et al. 
(1997, G97). Points from G97 have been recahbrated as described 
in the text. Frequency error bars on ACT and G97 points indicate 
bandwidth. Temperature error bars on G97 points do not include 
cahbration uncertainty. The green dashed hne is a spectral model 
from Weisstein &; Serabyn (1996, WS96), which shows proximity 
of absorption features to G97 points. 



Figure 6. Binned data and resulting best-fit two-component 
model of effective Saturn brightness vs. ring opening angle (B) 
relative to observer. The model is symmetric about B = by con- 
struction. As the absolute value of B increases, the increased radia- 
tion from the rings and decreased flux from the obscured disk leads 
to a local minimum in the total effective brightness at B ^ 9° (13°) 
at 149 GHz (219 GHz). The model is fltted independently for each 
frequency band. Error bars correspond to the error in the mean of 
observations contributing to each point, but flts include additional 
error due to calibration, which is covariant within each season and 
in some cases between seasons. Approximate mean dates of obser- 
vation are, from left to right: Nov/08, Dec/08, Jun/10, Apr/10, 
Dec/09, Dec/10, Dec/10. The sign convention is such that B < 
corresponds to negative values of the sub-Earth latitude. 



In Figure 7 we show the ACT disk temperature mea- 
surements along with high precision results from Wll 
below 100 GHz and from Goldin et al. (1997, hereafter 
G97) between 170 and 700 GHz. The G97 measurements 
have statistical errors at the 2% level, but are calibrated 
with reference to the same Mars model used by G&O. 
We have applied a rough recalibration factor to the G97 
measurements using the WMAP result (Weiland et al. 
2011) that the Ulich (1981) model for Mars temperature 
is 5% high at 94 GHz. Our recalibration factor varies 
linearly with the log of frequency from 0.95 at 90 GHz to 
1 at 857 GHz. With this calibration factor applied, the 
G97 measurement at 170 GHz and ACT measurements at 
149 and 219 GHz are consistent with a flat spectrum over 
this frequency range. It is difficult to comment further 
on the continuum spectrum of Saturn based on the ACT 
and G97 data, due to the proximity of the high frequency 
G97 bands to PH3 resonances (e.g., Weisstein & Serabyn 
1996). 

6. CONCLUSION 

We have described the reduction of ACT planet ob- 
servations for the purposes of measuring and monitoring 
the telescope beam and gain parameters. The core of 



the beam has been well characterized by Saturn obser- 
vations, and studies of point sources in the map have 
provided an empirical correction to account for pointing 
variance, alignment error, and slight changes in beam fo- 
cus. Observations of Uranus have been used to measure 
and account for changes in system gain due to atmo- 
spheric water vapor and focus variation. 

Based on hundreds of brightness measurements of 
Uranus, we have obtained precise measurements of the 
disk brightness temperature at 148 GHz and 218 GHz. 
In contrast to previous observations, our absolute 
calibration standard is the CMB blackbody. We have 
also obtained precise measurements of the Saturn disk 
temperature, in the context of a simple two-component 
model. 
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